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ABSTRACT 

The rate at which interstellar gas is converted into stars, and its dependence on environment, 
is one of the pillars on which our understanding of the visible Universe is build. We present a 
comparison of the surface density of young stars (S*) and dust surface density (Sjust) across 
NGC 346 (N66) in 115 independent pixels of 6x6pc^. We hnd a correlation between E* and 
Edust with a considerable scatter. A power law ht to the data yields a steep relation with an 
exponent of 2.6±0.2. We convert Edust to gas surface density (Egas) and E* to star formation 
rate (SFR) surface densities (Esfr), using simple assumptions for the gas-to-dust mass ratio 
and the duration of star formation. The derived total SFR (4±110“^ Mq yr“^) is consistent 
with SFR estimated from the H a emission integrated over the H a nebula. On small scales 
the EgFR derived using H a systematically underestimates the count-based Esfr, by up to a 
factor of 10. This is due to ionizing photons escaping the area, where the stars are counted. 
We find that individual 36 pc^ pixels fall systematically above integrated disc-galaxies in the 
Schmidt-Kennicutt diagram by on average a factor of ^7. The NGC 346 average SFR over a 
larger area (90 pc radius) lies closer to the relation but remains high by a factor of ~3. The 
fraction of the total mass (gas plus young stars) locked in young stars is systematically high 
(^10 per cent) within the central 15 pc and systematically lower outside (2 per cent), which 
we interpret as variations in star formation efficiency. The inner 15 pc is dominated by young 
stars belonging to a centrally condensed cluster, while the outer parts are dominated by a 
dispersed population. Therefore, the observed trend could reflect a change of star formation 
efficiency between clustered and non-clustered star-formation. 

Key words: ISM: individual: NGC 346 - ISM: Structure - stars: pre-main-sequence - galax¬ 
ies: Magellanic Clouds 
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1 INTRODUCTION 


It is generally accepted that the local star-formation rate (SFR) 
scales with the local gas density of the interstellar medium (ISM). 
This wisdom is captured in the Schmidt-Kennicutt relation (SK- 
relation, see for a review [Kennicutt & Evans|2012[ l which connects 
the observed gas-surface density (Egas) along various sight-lines 
with the observed surface density of SFR (Esfr). Even though 
there are multiple incarnations of the Schmidt-Kennicutt diagram, 
depending on which tracers are used to measure both Egas and 
Esfr, and whether to consider only the dense (molecular) gas or 
include the atomic phase, the emerging picture is clear: the higher 
the ISM column density, the higher the star-formation rate. A key 
feature of the SK-relation is its power law behaviour over a large 
dynamic range of gas column-densities and star-formation rates. 
The slope of the relation depends on whether one measures only the 
dense gas (exponent ~1) or also the atomic gas (exponent ~1.4). 


The physical interpretation of this interdependence, i.e. how 
to derive the efficiency with which interstellar gas is converted into 
stars, is much less self-evident, especially on small spatial scales 
(<1 kpc). This is due to the following: i) The ‘star-formation rate’ 
measures stars that have already formed, while the gas column den¬ 
sity relates to (the capacity for) future star formation. On small 
scales, with low numbers of stars, one can no longer expect to 
measure a time-averaged efficiency as is the case when averaging 
over a large number of star forming sites, ii) Dynamical simula¬ 
tions (e.g. [Kiseleva et al.||7998^ and observations (e.g. [Goodman] 
[2004 [ > show that young stars may be displaced, by several pc per 
Myr, from the location where they were bom, and thus removed 
from the reservoir of gas from which they formed Hi) The large 
scale averages include non-star-forming parts of the ISM and the 
relative proportion of this ‘inactive’ gas will change when zooming 
in on smaller scales, iv) Many of the used SFR-tracers (e.g. Ho? 
or far-ultraviolet (FUV) emission) are sensitive only to high-mass 
stars, which represent only a fraction of the stellar mass assembled 
in the star-formation event and to derive the total SFR one needs to 
extrapolate the mass-function to lower masses. Since the high-mass 
slope of the initial mass function (IMF) is very steep - very few of 
the formed stars are high-mass - the corrections are large and sen¬ 
sitive to the assumed shape of the IMF. Some or all of these effects 
are probably the reason why the studies of star formation rates in 
nearby individual molecular clouds i Evans^^etalJ|2009| [Wu et ak] 


|2010[[Heiderman et al.|201^|Gutermuth et al.|20i 1^ find very dif¬ 

ferent star formation rates at a given gas column density compared 
to global averages. 


The aim of this paper is to circumvent some of these compli¬ 
cations by studying the correlation between newly formed stars and 
the gas reservoir (Esfr, and Egas) in many small beams (~6pc) 
which cover an entire star forming complex. We investigate, thus, 
the small scale relation between young stars and gas reservoir while 
covering a larger region. We use the most direct way of determining 
Esfr, i.e. by counting the newly assembled stars. At the same time 
we employ the ‘dust method’ to determine the ISM column density, 
i.e. using the infrared (IR) to submillimetre (submm) dust contin¬ 
uum emission to derive dust surface densities. This method has the 
advantage of not being sensitive to the state of the gas (ionized, 
neutral or molecular) as long as dust and gas are well mixed. This 
insensitivity is important because we will use a catalogue of low- 
and intermediate-mass pre-main sequence (PMS) stars plus upper- 
main sequence (UMS) stars, i.e. O and B type stars, as the basis for 
our star counts. Due to the evolutionary stage of the PMS stars - a 
few Myr after the onset of star formation - the parental molecular 


complex may have been subject to significant photo-dissociation 
and using a molecular gas tracer (such as CO) may lead to an in¬ 
complete view of the total gas surface densities. 

The object of choice is the star-forming complex related to 
NGC 346, also often referred to as N66. NGC 346 is the most 
prominent star-formation region in the Small Magellanic Cloud 
(SMC). In order to assess the Esfr, an unbiased census of the 
young stellar population is required. This is available for NGC 346 
thanks to HST observations ( [Nota et ^[2006] [Gouliermis et al.j 
[2006] see Sect. O- We have access to sensitive IR to submm 
photometry from Spitzer (PI Gordon), Herschel (PI Meixner) and 
APEX/Laboca (PI Hony) in order to derive reliable dust surface 
densities. The proximity (60 kpc) of the SMC in combination with 
the size of the star forming complex (~100x 100 pc^) enables us 
to ‘zoom in’ (~6pc linear scale) and study local variations at the 
angular resolution of the available submm observations. 

Moreover, the entire complex encompasses a range of envi¬ 
ronments. In terms of the ISM, the region exhibits a wide range of 
surface brightness of hot dust (24 pm), polycyclic aromatic hydro¬ 
carbon (PAH; 8 pm) and molecular gas (CO(J=2-l), [Rubio et al.| 
|2000[ >. In terms of the stellar content, there are zones which are 
cluster dominated as well as regions where the young stars are more 
dispersed (Gouliermis, Hony & Klessen|2014[ l. This allows us to 
examine how the environment influences the relation between star 
formation and the available ISM reservoir. The main questions we 
address are: 


• How well do young stars and the densest ISM track each 
other? 

• Do star counts yield the same Esfr as indirect tracers? 

• Does the efficiency of star formation vary between clustered 
and more extended spatial distributions? 

• What is the effect of averaging over different scales on the 
derived parameters in the position in the SK-diagram? 


Additionally, we examine whether there may be an effect of the 
reduced metal abundance of the SMC (0.2 solar, |Haser et al.[1998[ ) 
on the star formation efficiency. 

The paper is organised as follows. Previous work on deriv¬ 
ing SFR from individual young sources is reviewed in Sect, 

In Sect. 1^ we present the catalogue of young stars and the IR to 
submm data-sets that we use for our analysis. We further describe 
the data treatment to derive independent surfac e den sity measure¬ 
ments in 20x20arcsec^ (6x6pc^) beams (Sect. 3.4 1 . Our method¬ 
ology to obtain gas surface densities, using SED fitting, is described 
in Sect. [4.1[ and [4.2[ The SER derived from our young star counts 
is compared to indirect SFR tracers in Sect. [4.3[ In Sect.J^we com¬ 
pare the derived Egas and Esfr- Possible biases that may affect the 
observed value of EsFR/Egas are discussed in Sect.]^ We conclude 
in Sect. [7] 


2 PREVIOUS WORK 


Most of the work done on comparing ISM column densities and 
star-formation rates based on individual objects has concentrated 
on nearby (within one kpc from the sun) Milky Way star-forming 
sites (Evans et al.||2009[ [Lada, Lombardi & Alves||2010[ |Heider-| 


[man et ^|2010[ [Wu et af]|2010[ Gutermuth et al.|[2011( . These 
studies have consistently shown that when one considers the im¬ 
mediate dense environment (on a few parcec scales) surrounding 
the newly formed stars, the derived SER is higher than predicted 
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Figure 1. False colour composite image (green: optical, red: infrared, blue H a) of NGC 346 based on HST imaging (PI Nota, |Nota et al.|2006 ). The contours 
of the, 2 arcsec resolution, density map (S*) of the young stars are shown in yellow. The contour levels correspond to 3, 10 and 20 stars pc“^. The maximum 
value of E* is 45 stars pc“^. 


by the galaxy averaged SK-relation (by up to a factor of 54) for a 
given ISM column density. 

[Evans et al.| ( |2009) ; |Lada, Lombardi & Alves] ( |2010| l present 
the total star-formation rates based on young stellar object (YSO) 
candidates, identified from their mid-IR excess using Spitzer and 
compare those to molecular cloud masses derived using extinction 
maps for ~10 nearby molecular clouds. These studies show that 
there is a large variety of star formation efficiency (SEE) among 
those clouds. [Lada, Lombardi & Alves| j2010^ show that some of 
the observed spread in SEE disappears when considering only the 
mass of gas above an extinction threshold of Ax > 0.8 mag. They 
interpret this as an indication for an density threshold below which 
star-formation does not occur. A similar conclusion is reached by 
[Heiderman et al.| ( |2010^ who study the Esfr versus Egas within 
those clouds using YSOs. The typical spatial scale over which they 
compare these quantities is from a few tenths of a parsec to a few 


parsec. Their data show a very sharp decrease in Esfr below Egas 
Mq pc“^ which is interpreted as resulting from a density 
threshold for star-formation. Note, that such a threshold, if it exists, 
cannot be a universal threshold. There are several observations of 
specific regions whose gas surface densities exceed the threshold, 
with much less star formation than would be expected (e.g.|Long-| 
[more et al.|2013[|Rathbome et al.|2014) . Moreover, the existence 
of a column density-threshold for star formation in a turbulent ISM 
is not supported by theoretical models ( [Clark & Glover|2014^ . 

A different picture emerges from the analysis presented by 
[Gutermuth et al.| ( |2011^ . These authors present scatter diagrams of 
surface density of from YSOs versus Egas (from extinction) for 
eight molecular clouds within 1 kpc from the sun. They find a steep 
power law correlation (with indices between 1.4 and 2.7, with an 
outlier at 3.8) with significant scatter and no indication for a density 
threshold. They argue that the data are consistent with a SLR which 
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Table 1. Values and symbols used throughout this paper. 


Quantity 

Symbol 

Value 

Comments/Refs 


SMC distance 

dsMC 

60ib3 kpc 

Harries, Hilditch & Howarth 20031 

Detected young stars 

-^star 

5150 

Gouliermis et al.|' 

2UU6^, see Sect. 3.1 

Total young stellar mass 

Mtot 

2.2-104 ]y[^ 

Sabbi et al. 

2D08 

, see Sect.|3.l| 

Mass per catalogue source 

Meat “ 

4.3 M© 

= Mtot -Y-.L, see Sect.|3.1| 

SF duration 

AtsFR 

5-10® yr 

Mokiem et al.' 2006|, see Sect.|3.l| 

Gas-to-dust mass ratio 

r-gd 

1250 

Sect. 4.2 




Derived Quantity 


Young star surface density 

s* * 

Stellar mass surf. dens. 

Em. 

SFR surf. dens. 

Esfr 

Dust surf. dens. 

Edust 

Gas surf. dens. 

Egas 

Stellar mass fraction 

fraCM- 


from star catalogue 

— ■ ^ITcat 
= Sm* 

from SED fitting 

“ ^dust ■ t*gd 

— (^M* + ^gas) ^ 


“This mass is not the mean mass of the HST detected sources but the mass each source represents after correcting for completeness. The mean mass of the 

young stars in the HST catalogue is ~2 Mq . 


is proportional to Egas- Note that the different conclusion reached 
by the latter authors is based, in part, on the way the surface density 
of YSOs is calculated in the low surface-density regime as they use 
the nth-nearest-neighbour centred on the detected YSOs. We fur¬ 
ther note that, given the steep dependence on Ssfr and the low 
number counts, distinguishing between a threshold and this power 
law is difficult and almost a matter of semantics. 

The work by [Gutermuth et al.| ( |2011^ is the Galactic study 
that can be compared most directly to the NGC 346 analysis we 
present here, because their Egas maps cover several tens of parsecs 
around the molecular clouds and each cloud harbours a significant 
number (>100) of YSOs. Compared to the work we present here, 
these previous studies have the advantage of proximity and as a 
result, they are able to examine the star-formation within individual 
molecular clouds on spatial scales of < 1 pc while the best angular 
resolution we can attain in the SMC in the ISM tracers is ~6 pc. 

Recent Spitzer surveys of the Large and Small Magellanic 
Clouds (LMC and SMC) made it possible to use resolved massive 
YSOs to determine current SFRs outside of the Galaxy. Unlike the 
consistently higher SFRs of small-size clouds reported in the afore¬ 
mentioned Galactic studies, [Chen et al^j2010[ > showed that current 
SFR determined from massive YSOs in individual giant molecu¬ 
lar clouds in H II complexes in the LMC are distributed above 
and below the SK relation, depending on the evolutionary states 
of the giant molecular cloud. The scatter in SFR diminishes and 
follow the SK relation when properties are averaged over the en¬ 
tire H II complexes of 150—200pc. This is in agreement with the 
findings of jKruijssen & Longmore| ( [2014^ , who attribute the scatter 
on smaller scales to incomplete statistical sampling of independent 
star-forming regions. Furthermore, the commonly used SFR trac¬ 
ers such as H a and 24 pm luminosities would under-estimate the 
amount of star formation in regions not actively producing 0-type 
stars, such as the LMC’s largest molecular ridge and the Magel¬ 
lanic Bridge. They show that using H q or 24 ^m as SFR tracers 
these points are outliers on the SK-relation while they follow the 
SK-relation when using SFR determined from YSO counts pnde-j 
jbetouw et al.|200^|Chen et al.|20l4T l. 

The work presented here differs from these previous stud¬ 
ies by the sources we use to trace the star formation (PMS-l-UMS 
stars versus YSOs) and the way we derive the ISM column den¬ 


sities (dust emission versus extinction). This has the following 
(dis-)advantages: 

• Rich stellar catalogue: //SThas provided accurate photometry 
for more than 5 000 stellar sources covering the pre-main-sequence 
( jCouliermis et al.|2006| >. 

• Measured mass function and age: ages and the mass function 
have been determined for the UMS-l-PMS stellar catalogue jSabbij 
jet al.|200'8| l. This means we can rely on the observed mass function, 

which covers a significant (0.8—40M©) part of the entire mass 
range, and the observed duration of the star-formation event to de¬ 
rive the SFR. The YSOs studies assume a typical age and mass for 
a YSO, both of which are uncertain (e.g. |Preibisch||2012^ , which 
could introduce a large systematic uncertainty. 

• Large coverage: the IR and submm observations that we use 
to measure dust column densities are very sensitive and most cover 
the entire SMC (Sect. [T^. As a result, we have access to an un¬ 
biased probe of the dust surface density out to large distances and 
low column densities. This is important to assess the behaviour of 
the SFR at low Egas, where significant disagreement exists in the 
literature (see above). 

However, since we are using stars that are no longer embedded, 
the sources in our catalogue are more evolved (<5 Myr, Sect. [13 
than YSOs (<2Myr) and have had more time to displace (or be 
displaced from) their parental gas reservoir. Moreover, since the 
whole region has to be more evolved for it to host so many PMS 
stars, the radiation from the hot, young stars has had more time to 
affect the surrounding medium. 


3 OBSERVATIONS 
3.1 Young stellar catalogue 

Imaging with the Hubble Space Telescope {HST) of star-forming 
regions in the Magellanic Clouds has provided access to excep¬ 
tionally rich samples of newly-bom stars down to the sub-solar 
regime over large areas of the sky. Our catalogue of young stars in 
NGC 346, one of the largest star-forming complexes in the SMC, 
is among the richest Magellanic Clouds stellar ensembles ever col¬ 
lected with HST. Three pointings (with significant offsets), centred 
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on NGC 346 were observed with the Wide-Field Channel (WFC) 
of ACS covering an area of about 5x5arcmin^ (~88x88pc^). 
The images were obtained within the FIST Program GO-10248 (PI 
Nota, |Nota et al.|[2006^ in the filters F555W and F814W (equiv¬ 
alent to standard V and I bands). Detailed descriptions of the data 
reduction and photometry are given in |Gouliermis et al^j2006^ . The 
photometry reveals more than 98 000 stars at various evolutionary 
stages, and its 50% completeness limit is msss ~ 27 jGouliermis| 
let al.|2006| l. 

The sample of young stellar populations in NGC 346 consists 
of the low-mass pre-main-sequence (PMS) stars, which are identi¬ 
fied from their faint-red positions on the colour-magnitude diagram 
(CMD). See |Gouliermis| ( (T012| ) for a review on low-mass PMS stars 
in the Magellanic Clouds. This sample is complemented by the 
UMS stars (selected with colours msss — msi 4 ^ 0.0 and mag¬ 
nitudes 12 ~ msss ^ 17), including more than 30 early-type stars 
(see also [Massey, Parker & Garmany|198^|Mokiem et al.|2006| >. 
The masses of the selected UMS stars are such that they would have 
evolved off the main sequence if they were older than the PMS pop¬ 
ulation. Thus, the catalogue we use in this analysis contains 5 150 
young (PMS and UMS) stars in the mass range between ~ 0.5 Mq 
and ~ 40 Mq . The location of the young stars within the complex is 
indicated in Fig.[T]with yellow contours. This image shows the very 
prominent H a emission associated with the central cluster and the 
presence of several secondary clusters of young stars to the North 
which cause little H a emission. 

It is hard to get accurate age and mass determinations of the 
individual PMS stars based on their position in the CMD jjeffries| 
|2012| l. An isochronal median age around ~ 3 Myr has been de¬ 
fined for the ensemble of PMS stars in NGC 346 on the basis of 
comparison with evolutionary models ( |Sabbi et al.|2008 | |Cignoni| 
|et al.|20rT| ). |Mokiem et al.|p006| l have spectroscopically studied 
the high-mass stars towards NGC 346. They find ages between 0 
and 5 Myr for the high-mass stars that belong to the star-formation 
complex. This age range is fully compatible with the position and 
spread of the low-mass PMS stars in the CMD. Since our objective 
is to measure the average SFR counting all stars that have formed 
during the star-formation event, we use 5 Myr as the star-formation 
duration (Afspu). Thus, AfsFR refers to the duration of the event 
not to the ages of the individual stars. Therefore, this approach is 
conceptually different from the YSO-counting studies (e.g. |Lad^ 
|Lombardi & Alves|2010[ |Heiderman et al.|2010| l in which ages of 
the individual objects have been assumed. 

We construct the entire mass function of the UMS-l-PMS cata¬ 
logue by extrapolating the mass function for masses above 0.8 Mq 
derived by |Sabbi et^T] ( |2008[ Fig 3 bottom panel) from the same 
observational material. Their mass function is completeness cor¬ 
rected for crowding and photometric sensitivity but not for inter¬ 
stellar extinction. Over the measured range the mass function has 
a slope of -1.430 We extrapolate to lower masses using a Kroupa 
mass function with a slope of -0.3 between 0.08 and 0.5 Mq and 
a slope of 0.7 between 0.01 and 0.08 Mq ( |Kroupa|2001| >. The to¬ 
tal stellar mass formed in the last 5 Myr, obtained by integrating 
this mass function and multiplying by the survey area (88 x 88 pc^) 
equals ~2.24' 10"* Mq. Since we have 5 150 stars in our catalogue, 
each source represents on average a mass of 4.3 Mq (Meat). The 


^ We use the formalism of the mass function - logio(Af(Alogio(m))) ver¬ 
sus logio(m) - as given by |ScaloH1986) . In this definition the Salpeter mass 
function has a slope of -1.35 and the Kroupa mass function has a slope of 
-1.3 above 0.5 Mq 


mass in directly detected stars (M* > 0.8Mq) is 1.2-10^ Mq, i.e. 
50 per cent of the total. The main parameters used in our analysis 
are summarised in Tab.[T] 

The total star formation rate, averaged over the last 5 Myr, of 
the complex calculated in this way is 4±1'10~^ Mq yr“^, where 
the uncertainty takes into account possible star formation outside of 
the field-of-view of HST (20 per cent) and the extrapolation of the 
observed mass function to lower masses (20 per cent) but not the 
systematic uncertainty due the assumed duration of the star forma- 
tion ev ent. This value is consistent with the results of ISimon et al.l 
(2007 I who find a SFR >3.2-10“® Mq yr“® over the last 1 Myr 
using Spitzer identified YSO candidates. 

For the rest of the analysis we assume that the mass function 
and the duration of star formation are invariant across the star form¬ 
ing complex. There are some indications of variations in the slope 
of the mass function in NGC 346 ( jSabbi et al.|2008||Cignoni et al.| 
|2011| |. The general trend is a steeping of the slope with increasing 
distance to the central cluster, i.e. a larger fraction of high-mass 
stars towards the central cluster, consistent with a certain degree of 
mass segregation in the cluster. However, such variations have little 
impact on our results which are based on star counts, while lumi¬ 
nosity based SFR tracers (e.g. H a) are very sensitive to the num¬ 
ber of high-mass stars and thus the high-mass slope of the mass 
function. As an illustration of this invariance consider the standard 
Kroupa mass function ( |Kroupa|2001| l with a high-mass power law 
slope of -1.3. When we vary this slope between -1.0 and -2.0 - the 
observed range of slopes in NGC 346 ( |Sabbi et al.|2008^ - the total 
mass per number of stars above 0.8 Mq (i.e. Meat) varies by only 
~20 per cent. For the same variations in slope the fraction of mass 
locked in the high-mass stars - those powering the H a emission 
(M* >8 Mq) - varies by 1500 per cent! 


3.2 Herschel and Spitzer data 

We use the mid-IR to submm photometric maps that have been 
obtained with Spitzer (SAGE-SMC :|Gordon et al.|2011) and Her¬ 
schel (HERITAGE: |Meixner et al.||2013^ . Within these two large 
programs the entire SMC has been observed at IRAC (3.6, 4.5 5.8 
and 8.0 ^m), MIPS (24, 70, and 160 pm), PACS (100 and 160/mi) 
and SPIRE (250, 350 and 500 pm) wavelengths. The sensitivity of 
these data sets is such that for the environment of the prominent star 
forming region NGC 346 the photometric uncertainty is dominated 
by the absolute flux calibration (see also Sect. |3.4^ . An overview of 
NGC 346 in some key ISM tracers is shown in Eig.0 In the follow¬ 
ing analysis we mainly use the IRAC, MIPS-24, MIPS-70, PACS- 
100, PACS-160 and SPIRE-250 data, because these sufficiently 
cover the spectral energy distribution (SED) to derive accurate dust 
surface densities at optimal angular resolution (~20arcsec). The 
SPIRE 350 pm and 500 pm data have poorer angular resolution. 
The only exception is the dust masses derived within large circular 
apertures (Tab.0 which uses the SPIRE-350 pm data to replace the 
Laboca because the latter is not observed over such a large area. 


3.3 Laboca 870 pm imaging data 

To constrain the emission by the coolest dust, we use Laboca 
data at 870/rm, instead of the SPIRE 350 data. This data set al¬ 
lows us to obtain measurements at ~20 arcsec resolution which is 
unattainable with HerschellSPlRE 350 and 500. Laboca is a multi¬ 
channel bolometer array for continuum observations at 870/rm in¬ 
stalled on the Atacama Pathfinder Experiment (APEX) telescope 
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Figure 2. False colour image showing an overview of the large environ¬ 
ment of NGC 346 in key ISM tracers. We show 8 /im emission indicative 
of the location of the most prominent PDRs (blue; IRAC4 emission minus 
the stellar continuum estimated from IRAC2), the warmer dust traced by 
PACS 100 (green) and the cooler dust (SPIRE 250, red). Each waveband 
is shown at its native resolution. The blue contours show the stellar density 
(see Fig. [TJ. The general observed trend is a decreasing dust temperature 
with distance to the most rich cluster (yellow to red colour gradient) and 
strong PDRs near (but not peaking on) the stellar concentrations. The re¬ 
gion where both dust surface density, including Laboca measurements, and 
stellar densities are well defined is indicated by the white line (see Sect. |3.4| 
for details). 



in Chile. It has a field-of-view of 11.4arcmin and a PSF full- 
width-at-half-max of 19.2±0.3 arcsec. Observations of the NGC 
346 region with Laboca were carried out in Decemher 2008 and 
April, May and September 2009 (Program ID: O-081.F-9329A- 
2008 PI: Hony). The mapping was done with a raster of spiral 
pattern in order to fully sample the final map (size of the region 
mapped: 25x25 arcmin^). Data are reduced with the BOlometer 
Array Analysis Software (BoA |Schuller et al.|[2010^ . The scans 
are reduced individually. They are calibrated using the observa¬ 
tions of planets (Mars, Neptune, Uranus, Venus) and secondary cal¬ 
ibrators (G10.62, G5.89, HLTAU, PKS0506.61, PMNJ0210-5101, 
PMNJ0303-6211, PMNJ0450-8100, N2071IR and VY-CMa). We 
determine zenith opacities using a linear combination of the opac¬ 
ity determined using sky-dips and that computed from the precip- 
itable water vapo^ We remove dead or noisy channels as well as 
stationary points or data taken at fast scanning velocity or above 
an acceleration threshold. We subtract the correlated noise, apply 
steps of median noise removal, baseline correction and despiking. 
The reduced scans are then combined into a final map. We isolate 
the pixels above a given signal-to-noise ratio to define a ‘source 
model’ that is subtracted from the time-ordered data. We rerun the 
reduction pipeline and add the source model to the time-ordered 
data to obtain a new map. The process is iterated until the result 
map converges (5 iterations were applied for NGC 346|^ The final 
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Figure 3. Laboca 870 /im map of NGC 346. The contours show the 36 cm 
radio continuum emission from SUMSS jMauch et al.|2003) . The 3 and 
6 cm radio continuum images (not shown) exhibit the same moiphology 
jPickel et al.|2010) . The possible corresponding contribution of free-free 
emission at 870 fim has been filtered-out by the median filtering applied in 
the data reduction process. The dashed circle indicates the approximate ex¬ 
tent of the H a nebula. We show a close-up of the brightest region in the 
inset, in which the 250 fzm dust emission contours, in perfect correspon¬ 
dence with the Laboca map, are shown. 


map is shown in Fig.[^ Within the region of interest there are three 
main emission features. There is an emission peak towards the cen¬ 
tral cluster, a direct indication that this sight-line is dust rich. We 
detect a brighter emission peak towards the north-west of the main 
cluster and a plume of 870 /rm emission which follows the string 
of stellar clustering towards the north-west. This plume coincides 
with the brightest CO(j=2-l) emission detected in the region [^bio| 
|et al.| ( [2000] >. Due to the median filtering that is applied to remove 
the atmospheric and telescope signal, any spatially very extended 
emission has been removed as well. This explains why we do not 
detect any free-free emission from the ionized gas which should be 
present in this region (see Fig.[^ 


3.4 Beam matched data 

In our analysis, we use maps with a resolution of 20 arcsec (6 pc 
at the distance of the SMC). All maps with intrinsic spatial res¬ 
olutions better than this resolution, including the stellar catalogue, 
have been convolved to the this effective spatial resolution. Laboca, 
MIPS 70 and SPIRE 250 maps have not been convolved since they 
already have an intrinsic spatial resolution very close to 20 arcsec. 
We use the convolution kernels available from Karl Gordons web 
pag^for convolving to the SPIRE 250 beams. Eor the stellar den¬ 
sity map - which has an exceedingly good intrinsic angular reso¬ 
lution we have directly convolved the 2 arcsec stellar density map 
with the SPIRE 250 /rm beam. After convolution, all maps have 
been re-projected on the same pixels scheme with pixels the size 


^ Tabulated sky opacities are provided at http://www. 
apex-telescope.org/bolometer/laboca/calibration 
More details on the techniques are provided in 


Galametz et al. 


2013 


http: // dirty .as.arizona.edu/~kgordon/mips/conv. 
psfs/conv_psfs.html 
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of 20 x20arcsec^ in order to have spatially independent measure¬ 
ments and to allow pixel-by-pixel SED extraction and comparison 
between the stellar densities and derived dust surface densities. We 
use the program swarp from AstroMatic]^ for calculating the re¬ 
projected maps, as it is in our experience the most accurate tool 
in conserving the mean surface brightness. This is particularly true 
when re-projectingto a pixel scheme with very different pixel sizes, 
which is our case[fj 

In order to create reliable random noise maps that represent 
well the noise in the original maps and the effects of convolution 
and re-projection we apply the following procedure: 

• The amplitude of the noise in each map has been estimated by 
measuring the median absolute deviation from the median (MAD). 
This value is a more robust measure of the noise in a map than 
the standard deviation in the presence of real source emission. The 
estimated standard deviation equals 1.4825*MAD. 

• We generate a pure noise map at the original resolution using 
the measured standard deviation. 

• We generate a projected noise map convolving and re¬ 
projecting this noise map in the exact same manner as we do the 
science map. 

• We measure the amplitude of the noise in the projected noise 
map as the standard deviation of the fitted Gaussian function to the 
histogram of pixel values. 

These output noise levels have been compared to the formal 
calibration uncertainties (expressed as percentages of the observed 
surface brightness levels). For all maps except the Laboca map, 
the absolute calibration is dominant throughout the region of inter¬ 
est. In the case of the Laboca map, we retain for each pixel in the 
projected noise map the larger value of the flux calibration uncer¬ 
tainty and the random noise value. These uncertainties should not 
be added (quadratically) because the Laboca flux calibration un¬ 
certainty reflects the scatter in the measured calibrator fluxes and is 
thus, in essence, also a random noise term. 

Pixels outside the area surveyed by HST and pixels where the 
Laboca measurement is less than 2 a are masked. There are 115 
valid, independent 20x20arcsec^ pixels (covering about ~4000 
pc^) in the following pixel-by-pixel analysis. Their outline is shown 
in Fig.|^ All maps at the other wavelengths have signal-to-noise ra¬ 
tios larger than 2 - in most cases much larger - in these pixels. 

4 DUST, GAS AND SFR SURFACE DENSITIES 
4.1 SED fitting 

In order to derive dust surface densities, we fit the IR to submm 
dust continuum SED for each independent pixel in the region. The 
dust emission fitting is performed with the |Galliano et al.| ( |2011a^ 
SED model. This model accounts for variations of dust excita¬ 
tion rates along the line-of-sight. Following the prescription of 
[Dale et ^ j2001^ , it assumes that the distribution of dust mass, 
Md, is parametrized by the starlight intensity, U, as a power law: 
dMd/dU oc U~°‘, between Umin and Umin + AU. The parameter 
U is the stellar flux received by the dust and is normalised to the 
solar neighbourhood value ( [Mathis, Mezger & Panagia|1983^ . The 

^ http://www.astromatic.net 

® The IDL/astrolib tool hastrom, often used for the same purpose, performs 
poorly in this case. 



Figure 4. Maps of the output parameters of the dust SED fitting. The region 
shown corresponds to the white contour in Fig. For reference we show 
the surface density map of young stars in contours (see Fig. We show 
the luminosity of the dust emission ranging from 150 to 1300 Lq (a), the 
average radiation field as experienced by the dust ranging between 6 and 
45 Go (b), dust surface density ranging from 0.05 to 0.3 Mq pc“^ (c) and 
the mass of PAHs relative to the total dust mass normalised to the Galactic 
value ranging from 0 to 0.3 (d). 


shape of the SED (in particular, both its peak wavelength and its 
width) constrain the parameters of this power law. 

For a given starlight intensity, U, the dust emissivity is derived 
from the jZubko, Dwek & Arendt|[2004| dust family: BARE-GR-S) 
model, designed to fit the diffuse Galactic emission. However, we 
replace graphite by amorphous carbon jZubko et al.|1996| AGAR). 
Indeed, [Galliano et al.| ( |2011a) showed that the standard Galactic 
grain properties applied to Spitzer and Herschel observations of the 
EMC lead to violation of elemental abundance constraints, while 
the modified composition (the ‘AC model’) gives realistic gas-to- 
dust mass ratios. In summary, we make the assumption that grains 
in the SMC will have similar optical properties as in the EMC. 

The choice of a set of grain emissivities is tainted by system¬ 
atic uncertainties. As demonstrated by [Galliano et al.| ( |2011b) , the 
derived dust mass change by a factor of 2 to 3 depending on the 
optical properties. However, it is also shown that the change is es¬ 
sentially a simple scaling of the dust mass, and do not significantly 
affect the trends. Besides, if the origin of the difference in emis¬ 
sivity between the Milky Way and the EMC is an evolutionary ef¬ 
fect, linked to metallicity (cf. [Planck Collaboration et al.|201l1 l, it 
is likely that using EMC emissivities puts us closer to the reality 
than simply assuming Galactic properties. 

The tit of the observed SED, with the previously described 
dust model, provides constraints on: the dust mass, Md \ the param¬ 
eters of the power law distribution of starlight intensities, Umin, 
AU and a, which we summarise by quoting the average starlight 
intensity, (U ); and the PAH mass fraction, normalised to the Galac¬ 
tic value, /pah. In order to estimate the final uncertainties on these 
derived dust parameters, we perform Monte-Carlo iterations, fol- 
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Table 2. Total dust mass, gas consumption time-scale and integrated SFE 
as a function of the area considered in the analysis. The first row corre¬ 
sponds to the area included in the pixel-by-pixel analysis. A radius of 52 pc 
is equivalent to the area covered by the HST observations. The 90 pc radius 
matches the prominent SPIRE 250 and 350 /rm waveband emission, while 
the largest radius probes the full extent of the H a nebula around NGC 346. 
The gas consumption time-scale (rgas) is defined as the Mgas/SFR and 
the integrated star-formation efficiency (SFEioMyr) calculated as the 
percentage of gas that will have been converted into stars after lOMyr, as¬ 
suming that the SFR remains constant. We use lOMyr as the limiting age, 
when supemovae might disrupt the ISM and halt the star formation. 


area/radius 

-^dust 

[Mq] 

'^gas 

[10® yr] 

SFE^OMyr 

[%] 

mask 

498 ± 18% 

1.7 

5.8 

52 pc 

900 ± 19% 

2.5 

4.0 

90 pc 

1616 ± 19% 

4.5 

2.2 

145 pc 

2741 ± 22% 

7.6 

1.3 


lowing the method described in |Galliano et al.| ( [201 la^ , taking into 
account the correlation between calibration uncertainties. The me¬ 
dian uncertainties on the derived parameters (a^lx) are 6 per cent 
for Lduat, and 19 per cent for Edust and (JJ). The median standard 
deviation on /pah is 0.02. 

Maps of the best fit output parameters are shown in Fig. 
One of the key results from the SED fitting which is evident from 
these maps is the clear distinction between the parameters which 
trace the heating of the dust (panels a and b) and the surface den¬ 
sity of dust (panel c). Clearly the dust luminosity in the very active 
star forming region NGC 346 (panel a) traces closely the radiation 
field. This is in stark contrast to other (more usual) observations 
of more diffuse regions where the incident radiation field (and dust 
temperature) does not vary much and Ldust traces the dust surface 
density. Panel d shows that the PAH abundance close to the main 
cluster is so low that their emission is undetectable in the broad¬ 
band photometry. This low abundance of PAH is in perfect agree¬ 
ment with the findings of |Sauvage, Vigroux & Thuanj ( |1990t who 
deduced an increased differential depletion of PAHs towards NGC 
346 based on IRAS photometry. Similar paucities of PAHs towards 
other high-mass star forming region in the Magellanic Clouds have 
been reported (e.g. |Sandstrom et al.|2010[|Hony et al.|2010| l. 


4.2 Conversion of dust surface density to gas surface density 


We assume that gas and dust are well mixed along each line-of- 
sight in which case the gas surface density (Egas) is proportional 
to the dust surface density (Edust) and the proportionality factor is 
the gas-to-dust mass ratio (rgd). The dust masses derived by SED 
fitting methods vary significantly and systematically (e.g. |Galliano| 
|et al.|2011b[ [Gordon et al.|2014^ depending on the applied method 
and assumed dust optical properties. It is therefore important to 
use a value of rgd consistent with the used SED fitting method¬ 
ology, i.e. the value of rgd that gives the correct gas mass given 
the derived dust mass. We have calibrated the rgd using the avail¬ 
able gas tracers on a large scale. Within a SOOarcsec radius cen¬ 
tred on NGC 346 we find an atomic gas mass of 2.7-10® Mq using 
the Parkes-l-ATCA HI map ( Kim et aL||2003| > and a spin temper¬ 
ature of 60K ( [Bernard et al.||2008^ . The molecular gas mass de¬ 
termined from CO(I=1-0) within the same aperture is 7.4-10"^ Mq 
when using a CO-to-H 2 column density conversion factor (Xco) 


of 10 


[cm ^ (Kkms ^) i Bolatto, Wolfire & Leroy|2013 i. If 


we use an Xco of 10^^ [cm“^ (Kkms“^)“^] we find a molecu¬ 
lar gas mass of 7.5-10® Mq. Such an elevated Xco may be more 
appropriate given the harsh radiation field that is prevalent in the 
region, which may induce enhanced photo-dissociation of the CO 
molecules. Thus, we derive a total gas mass in the range of 2.7-10® 
to 3.4-10® Mq. The dust mass derived using the dust continuum 
photometry in this aperture and applying the same SED fitting rou¬ 
tine we use for the pixel-by-pixel analysis yields 2.7-10® Mq. The 
derived rgd ranges from ~ 1000 to ~ 1250. In the following we use a 
fiducial value for Tgd of 1250, which is much higher than the Galac¬ 
tic value (100, |Draine & Li|2007) due to the lower metallicity of the 
SMC. Our value for rgd of 1250 is lower than the total SMC inte¬ 
grated value of 1740 from jCordon et aL| ( |2014| (. The global value 
includes a significant amount of atomic gas at low column density 
in the outskirts of the SMC which is not detected in dust emission. 
This drives up the global rgd. 

Given that the derivation of rgd is susceptible to significant 
systematic uncertainties in Xco for the gas surface densities and 
dust optical properties for the dust surface densities, we will dis¬ 
cuss the implications of a factor of two variations in the absolute 
value of rgd on our results. Note that, part of the purpose of this 
paper is to determine how well the ISM and stars follow each other 
on small spatial scales. This part of the analysis does not depend 
on the chosen absolute value for rgd. The comparison of SER and 
star-formation efficiencies with values derived for other regions is 
however sensitive to the chosen value and we will discuss the influ¬ 
ence of our choice on the results we find. 

We further assume rgd to be invariant for the entire region. In 
principle, it is possible that rgd could vary significantly, in partic¬ 
ular in the vicinity of high-mass star-formation sites, where dust 
destruction can become important due to supemovae explosions. 
[Reid et ar]j2006^ have studied the supemovae remnants within the 
large environment of NGC 346. These authors show that the region 
we study here appears largely unaffected by the impact of super¬ 
novae. The closest remnant they identified (SNR B0057-724, see 
also Eig. 1^ lies to the North-East of the main cluster, outside of 
the region of interest. The observational data argue against strong 
variations of rgd because they show that the ratio of dust heating 
compared to gas heating is very constant throughout the whole re¬ 
gion (see Sect.|6.2^. 


4.3 Star formation rate surface densities 


The star-formation-rate map has been derived using the observed 
number of young stars along the line-of-sight in each pixel of the 
beam-matched stellar density map. It is calculated as (see also Ta- 
ble[T](: 


[Mgyr-i] 


( 1 ) 


AfsFR 

The star-formation surface density map (Esfr) is obtained by di¬ 
viding this map by the area of a pixel in kpc®. 

In Eig.|^ we compare the derived Esfr with selected, indirect 
star-formation rate tracers. Note that we only compare with those 
tracers (Ehq- calibration from Haoeta^ 201 1^ for a Kroupa IME; 
E24 gm- calibration from jRieke et al. 1 2009 1 and Etir- calibration 
from jMurphy et al.|201 1^ which are appropriate for the very recent 
epoch of star-formation (<5 Myr) and we do not consider EUV be¬ 
cause its calibration assumes a constant SFR over >10—100 Myr. 
As can be seen, the Esfr from star counts is generally higher than 
those obtained from Ha or dust emission. Moreover, the differ- 
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Figure 5. Comparison of derived pixel-by-pixel SgpR using different 
methods. The x-axis shows the SsFR using the young-star surface den¬ 
sity. In blue we show Sspp calculated following |Hao et ^ j2011[ for 
the Kroupa IMF) using the continuum subtracted H a map form MCELS 
jSmith, Leiton & Pizarro|2000^ that was kindly provided by S. Points. The 
red crosses are Sspp estimates using the total infrared (TIR) (3—1100 /rm) 
surface brightness following |Hao et ar]j2011) and the black symbols show 
the estimated Sgpp based on Sphzer/MlPS 24 pm jRieke et al.|2009|. 


ences are largest where the stellar density (i.e. Espr) is highest. 
Towards the densest cluster the difference between counting stars 
and using Ha is a factor of 10, while the dust tracers yield a local 
Espr lower by a factor of 50 . The dust tracers yield much lower 
values (by a factor of ~ 5 ) than both H a due to the reduced dust 
abundance in the SMC. 

The origin of the weak dust emission in NGC 346 is not the 
same as the weak emission in the Gould-belt star-formation regions 
reported by | Vutisalchavakul & Evans|j2013^ . These authors have 
compared the SFR one would derive from Spitzer/MlPS 24 pm 
measurements or TIR with SFR derived from counting YSOs. They 
find that the IR tracers underpredict the true star formation rate 
by two orders of magnitude, because these regions lack high-mass 
stars due to incomplete sampling of the initial mass function. NGC 
346 is a very rich environment and does not have such a lack of 
high-mass stars. 

The difference between young star counts and the H a derived 
values is due to the fact that the individual massive young stars 
in NGC 346 produce H a (and dust) emission over a much larger 
surface than the size of the pixel in which they are counted, signif¬ 
icantly larger than the full HST field-of-view. Note, that this effect 
will cause the tight relation between SFR (derived from H a) and 
the gas column density to break down on small (,$100 pc) scales, 
similar to what is observed in M33 ( [Onodera et al.|[T010t . This 
break-down occurs even in the absence of drifting of the stars from 
their birth environment or spread induced by evolutionary effects, 
simply because the local H a (or dust) emission no longer measures 
the locally produced UV-radiation field. Figure|^demonstrates how 
large this distance is. We show the fraction of the SFR which is 
recovered when integrating the extinction corrected H a emission 
(based on the combination of H a and 24 pm emission) in apertures 
of increasing radius. Within the HST field-of-view we recover 50 


radius [pc] 



radius ["] 

Figure 6. Derived total SFR for the NGC 346 complex based on the com¬ 
bined emission of H q and 24 pm as function of the radius (centred on 
14.776°;-72.1735°) over which the flux is integrated. We follow the pre¬ 
scription from |Calzetti et alT]p007| Eqn. 5). The SFR is expressed in terms 
of the fraction of the total SFR derived from star counts. The H a and 24 pm 
emission are much more extended than the distribution of the young stars. 
As can be seen, the SFR derived within the radius of the HST coverage 
(-^50 pc) is about half of the total SFR. Only when we integrate over the 
complete nebula (>^500 arcsec), where the H a emission reaches the diffuse 
background level do we derive similar SFR values. 


per cent. When we integrate out to a radius of 500 arcsec (145 pc), 
where the H a emission from the NGC 346 complex becomes indis¬ 
tinguishable from the general SMC emission, we recover ~90 per 
cent. In this calculation we assume that outside the field-of-view 
of the HST observations, the star formation is insignificant. This is 
corroborated by the results of |Sewilo et al.| j2013[ l who identified 
YSO candidates based on IR photometry of the entire SMC. Of the 
sources from their highly-reliable YSO-candidate list which are as¬ 
sociated with NGC 346, less than 20 per cent (4/23) are located 
outside of the field-of-view of HST. A very similar result is found 
when using the spectroscopically confirmed early-type stars from 
[Bonanos et al.| ( |2010^ . Of the 43 O-type stars from this catalogue 
within a 100 pc radius around NGC 346 only ~15 per cent (6) are 
not within the field-of-view of HST. 

We conclude that the total SFR as measured from the (pre¬ 
dominantly low-mass) young star counts is consistent with the SFR 
measured from H a emission to within ~10 per cent. However, the 
absorption of the hard radiation from the hot stars and subsequent 
re-emission as either recombination lines or dust continuum occurs 
far - at distances > 100 pc, in the case of NGC 346- from the ac¬ 
tual location of these stars. See also |Relano et al.| ( |2012t andlLil 
|et al.| ( |2013^ for excellent discussions on how radiation escaping 
the local environment affects H a and dust emission as SFR trac¬ 
ers on small spatial scales. As a consequence, these indirect tracers 
which depend on absorption and re-emission, cannot be used for 
our present purpose, which is to relate the stellar densities with the 
ISM on smaller scales within the star forming complex. Note that 
most extra-galactic studies that use H a to derived SFR are probing 
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Figure 7. Scatter plot of surface density of young stars versus surface den¬ 
sity of dust per pixel of 20x20arcsec^ (6x6pc^). The dotted line shows 
the ‘standard’ relation from|Kennicutt & EvansH2012| Fig. 1 la). The values 
from lHeiderman et al.|^2010[ , based on YSO counts and ISM column densi¬ 
ties in nearby individual molecular clumps, are also shown. The values for 
NGC 346 (black pluses) lie systematically above the ‘standard’ line. The 
thick blue stars indicate where NGC 346 is located when we average over 
all pixels in the masked area, and subsequently over apertures with radii of 
52, 90 and 145 pc (from top-right to bottom left). These radii correspond to 
the approximate size of the HST surveyed area, the prominent far-IR emis¬ 
sion nebula and the size of the H a nebula, respectively. The power law 
exponent of the effect of averaging, i.e the slope of the an'ow in the figure, 
is 2.9. 


NGC 346 - individual pixels 
NGC 346 - averages (see caption) 
Gould belt - Heiderman et al. (2010) 
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much larger regions (^kpc) than the scale on which the indirect 
tracers are reconciled with the direct measurements (~ 150 pc) and 
are thus unaffected by this effect. 


5 COMPARISON OF GAS AND STARS 

In the following section, we compare the derived SFR and ISM 
column densities. We first concentrate on the global behaviour and 
how NGC 346 compares to other measurements before looking in 
more detail at the tendencies within the complex. 


which is currently, actively forming stars, while the galaxy aver¬ 
aged values include also regions which are not presently forming 
stars. Compared to the study from [Heiderman et al.| j2010^ , our 
approach has the advantage of using the dust in emission to mea¬ 
sure the column densities and that we have access to dust emission 
covering a large area surrounding the star forming complex. This 
allows us to calculate averaged values over relevant spatial scales, 
taking into account also parts of the ISM which are less vigorously 
forming stars. 

The thick blue stars correspond to averages over the entire 
masked area, the entire HST field-of-view, the prominent SPIRE 
250 ^m emission and the Ha nebula (top-right to bottom left, 
respectively, see also Tab. |^. The mean value lies progressively 
closer to the SK-relation as we increase the area over which we 
calculate the average. It should be noted that the same shifting to 
bottom left must occur also to the values from lHeiderman et al.l 
( [2010^ if they were to integrate to larger distance from their young 
stars (to lower column densities), as the conditions that prevail in 
their low-mass star forming regions are probably more represen¬ 
tative for the spiral galaxies at low Esfr rather than those on the 
dividing line between the most active spirals and star burst galaxies. 
The third star from the top in Fig.|^- highlighted by a circle - indi¬ 
cates the average over the area where significant cold dust emission 
is detected in excess to the general SMC emission. This is probably 
the best estimate of the average SFR and ISM column density for 
the complex in its entirety. This value is above the SK-relation by 
about a factor of three albeit roughly compatible with the average 
relation when taking into account its scatter of 0.3dex. This factor is 
consistent with the prediction from jKruijssen & Longmor^j2014| l 
for the size of the region ~300 pc diameter over which we do the 
comparison. We also indicate in Fig.|^by how much the NGC 346 
values would be displaced (left-to-right) relative to the data from 
[Heiderman et al.|(2010| > and the SK line as a result of a factor of two 
variations in rgd. If rgd were 2500 as opposed to 1250, the average 
values of NGC 346 would fall on the SK line. Such a high value for 
Tgd would imply that half of the total gas (more than 80 per cent of 
the molecular gas) within a 500 arcsec (145 pc) radius around NGC 
346 is molecular without corresponding CO emission, as a result of 
photo-dissociation of the CO molecule. Such an large reservoir of 
hidden molecular gas, not related to CO, seems difficult to reconcile 

[& Maloney[ ( [201l} show that the distribution of the [C II] 158 ^m 
line emission closely resembles the CO distribution. 


with the emission of C^ which traces the photo-dissociation. 


Israel 


5.1 Global comparison 

We show in Fig.j^the stellar surface density versus the dust surface 
density for each pixel. These are the basic observational quanti¬ 
ties which do not depend on assumptions about Meat, AIsfr or 
rgd. Using the values in Table^we convert these into Esfr and 
Egaa (right and top axes). The individual pixels within NGC 346 lie 
above the fiducial SK-relation for Esfr versus Ehi+H2 (|Kennicutt| 
[& Evans|2012| Fig. 1 la) by on average a factor of 7. For reference, 
we also show the points for individual clumps in nearby Milky Way 
star forming regions from [Heiderman et al.[j2010l l, which also lie 
consistently above the SK-relation. These authors use SFR from 
low-mass YSOs and column densities from the dense gas measured 
from extinction maps on a typical scale of a few parsec. The fact 
that all these small-scale measurements lie above the galaxy aver¬ 
aged values is expected ( [Kruijssen & Longmore[[2()14[ ). The off¬ 
set is due to the selection biased introduced by studying a region 


5.2 The physical origin of variations in stellar mass fraction 

On the spatial scale that we probe here (~6x6pc^) we measure 
significant scatter in the relation between Eg^g and Em* (Fig.|7]l. 
Such scatter on small scales is expected, due to the stochastic na¬ 
ture of the local star-formation process, as well as relative motions 
between the formed stars and the reservoir from which they formed 
( [Kru ijssen & Lon gmore|20l^ . The observed scatter is more easily 
assessed in Fig. where we concentrate on the values derived for 
NGC 346. This figure also shows the derived power law for the en¬ 
tire ensemble of sight-lines. We use the IDE routine MPFITEX^F 
([Will iams, Bureau & Cappellari|20ldl l which uses MPFIT [Marl^ 
[wardt[ ( [2009| l to determine a straight-line fit to the data (in log-log 
space) with uncertainties in both variables. The best fit value for 
the exponent is 2.6±0.2, i.e. the relation is much steeper that the 

^ http://purl.org/mike/mpfitexy 
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Figure 8. Scatter plot of surface density of young stars versus surface den¬ 
sity of dust per pixel of 20x20arcsec^ (6x6pc^) zoomed in on the range 
occupied by NGC 346. The blue points have a strong warm dust component 
in their SED (S 24 ^m/S 250 /rm [F^/Fi^] > 0.3). Clearly, the lines-of-sight 
selected using this colour criterion occupy a separate region of the diagram. 
The long and short dashed lines are the straight-line fits in logarithmic space 
to the data with or without the warm pixels, respectively. The red areas cor¬ 
responds to the 2(7 uncertainties of the best-fitting slope, while keeping the 
normalisation fixed. 

canonical SK-relation (exponent: ~1.4). This very steep increase of 
SsFR with Egas is similar to the large exponents found by|Guter-| 
|muth et 5I| ( |201lt . 

The observed variations in SFR at higher dust/gas surface den¬ 
sities are clearly correlated with other observables, in particular the 
dust colour-temperature (here E 24 ;jm/E 25 o^m [Jy/Jy]). Perusal of 
the data shows that those pixels with E 24 ;jm/E 25 o run >0.3, i.e. 
warm dust, are towards the rich stellar cluster and all exhibit high 
values in this diagram. The position of these warm sight-lines is 
also indicated in Fig. We further show the power law fit to the 
data excluding the lines-of-sight with this very warm dust (expo¬ 
nent 1.6±0.1). The derived exponent of the power law varies signif¬ 
icantly depending on whether we include these sight-lines, which 
is an indication that the lines-of-sight towards the main cluster oc¬ 
cupy a particular part of this diagram. 

It is interesting to note that the rms scatter around the best¬ 
fitting power law is ±0.4dex. This is similar to the rms observed 
(±0.3dex) in the SK-relation of integrated ‘normal’ spiral galax¬ 
ies jKennicutt & Evans|2012| ). Perhaps, some of the differences in 
physical conditions that underpin the variations among spiral galax¬ 
ies are encountered within the NGC 346 star forming complex. If 
this is indeed the case, then studies such as the one we present here, 
offer a rare opportunity to zoom in to specific environments and de¬ 
termine the parameters that drive such variations. 

In the following we discuss the variations in terms of the stel¬ 
lar mass fraction of young stars (/raC]y[,^ = Em* (Em* +Egaa)~^) 
to facilitate comparison with star cluster formation models, where 
Em* refers to the mass of the PMS-(-UMS stars. Moreover, we pre¬ 
fer this terminology over star formation efficiency, which concerns 
the time averaged rate with which a large quantity of interstel¬ 
lar matter (entire complexes) is converted into stars. On the small 


Right ascension [J2000] 

Figure 9. Map of the NGC 346 which shows the derived stellar mass frac¬ 
tion. The size of each point is dependent on /racM* and its colour depends 
on E 24 line The mass fraction is indicated in each point. It is clear from this 
diagram that the stellar mass fraction peaks on the stellar cluster where the 
dust is heated to higher temperatures due to the intense radiation field. The 
warm pixels identified in Fig.[^are encircled. 

scales we probe here, we measure the stars that have formed earlier 
and the remaining ISM column density. It is not evident that these 
stars have formed from the ISM we detect now along the line-of- 
sight nor that they will still be in this location, if and when the next 
generation of stars is bom. 

Across most of the surveyed area (66 per cent of the pixels) 
/racM* is below 2 per cent, with the most common value being 
1 per cenj^ see also Fig. However, there is a significant tail to 
much higher values with a maximum of 15 per cent. These very 
high values correspond to very high stellar surface densities, not to 
low ISM column densities (see Fig.|^. 

The observed variations in fracy^^ are certainly not random, 
but show a definite dependence on other measurable quantities. 
In particular, the highest values are found near the main cluster 
(Fig.|g. We also show the trend of stellar mass fraction as a func¬ 
tion of distance to the main cluster in Fig. Even though there 
are clear deviations from spherical symmetry in the region, we ob¬ 
serve a monotonic decrease in fracy^^ with distance which can be 
approximated by a power law with exponent -0.7. 

From the observations and the SED fitting, we have at our dis¬ 
posal a number of key probes that typify the conditions along the 
line-of-sight of each pixel. These range from quantities that mea¬ 
sure directly the stellar densities and the resulting intense radiation 
field, to those that depend on both the local radiation field and the 
ISM column density and state. In these same ordering, these probes 
are: E*, E24 ^m, the average radiation field (((/)) as experienced by 
the dust grainQ Egas, the molecular gas fraction and the molecu¬ 
lar gas surface density. Perusal of the data shows that fracy^^ cor¬ 
relates well with stellar densities and hot dust emission at 24 pm 


® 1 per cent corresponds to a gas consumption time-scale of 5T0® yr, given 
the duration of the star-formation event of 5 Myr 

® {U) is the average field weighted by dust mass and as a result does not 
measure the ambient radiation field if the line-of-sight extinction is signifi¬ 
cant. 
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distance [pc] 

Figure 10. Trend of stellar mass fraction as a function of distance from the 
main cluster (ra=14.755179°, dec=-72.171201° [J2000]). Individual pixel 
values are shown by red pluses. Annular-averaged values (in steps of 7 pc) 
are shown with black diamonds. The error bars shown are \<j uncertain¬ 
ties on the mean. The behaviour can be approximated by fracy^^ [%] = 
20%-distance [pc]^“^'^^ between and 55 pc (dashed line). 


(Fig. [m a,b). The correlation with the ISM column density is much 
less tight (panel c), which is in part due to the fact that the ampli¬ 
tude of variations in Egas are much smaller that those is E*, i.e. the 
gas map is much flatter than the star map. When we use Eco/Sgaa 
as a proxy for the molecular gas fraction (panel d), there is no obvi¬ 
ous correlation between the stellar mass fraction and the molecular 
gas fraction. Specifically, the highest and the lowest stellar mass 
fractions coincide with intermediate values of the molecular gas 
fraction. The same is true also for the molecular surface density 
(panel e). 

The basic conclusion we can draw is that the highest stellar 
mass fractions are found where the highest concentrations of stars 
are found, towards the central cluster. There are many stars where 
there are many stars. We find no obvious parameters of the cur¬ 
rent ISM towards the cluster that explain why the cluster formed 
there. It is neither the densest nor the most molecular ISM within 
the complex. Probably the strong radiation over the past few Myr 
has already erased the memory of the conditions that led to the for¬ 
mation of the cluster. 

It should be noted that the central location of the main clus¬ 
ter may be a clue to why it formed there. It truely is central to the 
complex in most tracers. In the ISM column density tracers, in par¬ 
ticular, the main cluster location appears as the intersecting ‘hub’ 
of three peripheral dense ‘arms’ (Fig.[2|and|12|(. If the main clus¬ 
ter is indeed located at the meeting point of several dense arms, it 
possible that the star formation process was/is being fed by these 
structures. Whether this is the case cannot be assessed without kine¬ 
matic measurements of the dense gas. 


6 POSSIBLE BIASES 

We find that the stellar mass fraction is highest towards the rich cen¬ 
tral cluster. There are several effects and biases that could influence 
this interpretation. 


6.1 Completeness 

The derived relation between the newly formed stars and the ISM 
depends on the completeness of the HST stellar catalogue. The av¬ 
erage completeness for stars which are too faint and crowding has 
been accounted for by integrating the completeness corrected mass 
function ( |Sabbi et al.|2008[ see also Sect. [^. If there are remain¬ 
ing incompleteness issues in the catalogue, the fracj^^ would be 
even higher that we derive now. The observed variations in fracy^^ 
could be affected, if there are large spatial variations in complete¬ 
ness. 

Crowding: The regions which could be affected by crowding more 
than the average are the regions of the highest stellar density. How¬ 
ever, these regions already have the highest mass fraction of stars, 
and an additional correction would place these regions farther from 
the mean rather than closer. 

Extinction: NGC 346 exhibits some sight-lines with significant 
amounts of dust (Fig. [m which could hide embedded clusters of 
young stars. This could affect, in particular, the regions labelled #3 
and #4. Both regions have a lower than average stellar mass fraction 
and therefore additional embedded young stars would bring them 
closer to the average. However, the data argue against a significant 
population of embedded sources. First, because the peaks in the 
Edust map do correspond to local peaks in the HST detected E* 
map. Secondly, because there are no corresponding peaks in the 
MIPS 24 ^m maps at these locations, while an embedded cluster 
would produce 24 fim emission. This emission would be detectable 
unimpeded by line-of-sight extinction. Finally, these regions have 
low observed dust temperatures, while embedded heating sources 
would increase their temperature. 


6.2 Variations in gas-to-dust mass-ratio 


If there are significant and systematic variation in rgd, the conver¬ 
sion from Edust to Egas would be affected. Specifically, if the dust 
in the central part has been mostly destroyed, we would grossly un¬ 
derestimate the Egas there, which would bring these points closer 
to the general trend (to the right in Fig. [^. The SED modelling 
shows that these regions are extremely depleted in PAHs. Is it plau¬ 
sible that the dust has also been destroyed? The increase required 
in the assumed value of rgd for the central region of NGC 346 in 
order to make the variations in fracy^^ disappear is a factor of 10 
(rgd equal to 12 500 instead of 1250). 

An argument against significant dust destruction (or more gen¬ 
erally significant variations in rgd) is the very similar behaviour of 
H a and TIR (and 24 fcm) emission. The rms of Etir/Ehci of the 
115 pixels across the whole region is only 20 per cent and its be¬ 
haviour is not spatially correlated with the stellar density or the 
radi ation field density. This similar behaviour is also reflected in 
Fig. f ° For the relatively low extinction conditions in NGC 346, 


The scatter in Sxir/Shq of 20 per cent is much smaller than the scatter 
- about the mean trend in Fig.[^- in the Stir and Siia, individually. This 
indicates that an above (below) average value of Sjia corresponds to an 
above (below) average value in Stir. 
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Figure 11. Scatter plots of the stellar mass fraction as a function of a number of key tracers. From left to right the tracer measures increasingly the ISM density 
and the presence of dense clumps. From right to left the sensitivity to the local stellar density increases. E* is the surface density of the PMS+UMS stars. All 
panels span an identical dynamic range (Xmax/atmin = S/max/j/min = 1 000). 


Ehq measures the chance for the starlight to he absorbed locally by 
the gas and Etir measures the same for the dust, the very similar 
behaviour (on scales of 6x6pc^) of these tracers suggest a simi¬ 
lar proportion of dust to gas across the region we studied. In case 
of very reduced dust abundances towards the central region, one 
would expect a limb-brightened shell in the dust emission, contrary 
to what is observed (Fig.|^). 

A relatively constant Tgd is further supported by the analysis 
of |Lebouteiller et al^ ( |2008^ who measured the gas phase abun¬ 
dances along 13 different sight-lines across the NGC 346 complex. 
They find only minor variations in the gas phase abundances. Most 
importantly they find only 0.1 dex variations in the S/H ratio while 
sulphur can be depleted into dust grains and strongly varying dust 
abundances would lead to variations of the abundance of sulphur 
remaining in the gas phase. We conclude that variations in Tgd that 
could explain the variations in frac^^ are improbable. 


7 DISCUSSION 
7.1 Effects of evolution 

Our observations provide a snapshot of NGC 346 at one instance, 
roughly 3—5 Myr after the onset of the main star-formation event. 
During the star-formation epoch the amount of gas transformed 
into stars must evolve. Is the observed scatter in fracy^^ caused 
by studying only this snapshot? In other words, is it possible that 
after the star formation has ceased, the fracy^_^ will be more homo¬ 
geneous? Within the complex, different regimes can be discerned, 
as is illustrated in Fig.|12[ 

• Region #1: Main cluster with high dust surface densities and 
high number of stars. The region is detected in CO with moderate 
brightness. 

• Region #2: Secondary cluster with fewer stars and less dust. 
Appears to be a ‘scaled-down’ version of the main cluster with 
bright H a and MIPS 24 /rm emission and a moderately high 
/’’“Cm*. 

• Region #3: This blob to the north-west of the complex has 
the highest dust surface density and low measured fracy^^. It is 
detected in CO but the estimated molecular gas fraction is low. 

• Region #4: The molecular ‘spur’ to the north-east has a small 
stellar mass fraction, high dust surface an no associated H a emis¬ 
sion. 


It is unlikely that these regions will all produce the same stellar 
mass fraction in the end. It can be argued that region #4, with its 
low current fracyi^ but large reservoir of molecular gas, has the 
potential to still form a significant number of stars and thus boost 
its final mass fraction. However, the conditions in region #1—#3 
do not seem to favour prominent continued star formation. Their 
ISM seems lacerated with ionized gas, as testified by the prominent 
H a emission and based on the CO map the molecular gas fraction 
is also low. We assess that, most likely, the observed differences 
between regions #1—#3 will persist. 


7.2 Gas removal and stellar motions 


We find a clear trend in stellar mass fraction from the central parts 
of NGC 346 to its outskirts. Within ~15 pc /rac,^j^ ~10 per cent 
while between 15-50 pc it is ~2 per cent. These zones are dom¬ 
inated by the cluster and the dispersed young stellar population, 
respectively as shown with the auto-correlation analysis of the stel¬ 
lar distribution ( [Gouliermis, Hony & Klessen|2014^ . The observed 
decrease of fracyi^ with distance from the cluster is a signpost 
of lower SFE on the periphery (dispersed SF) than in the densest 
central parts (clustered SF). However, there are two caveats to con¬ 
sider, both related to the relative displacement of the stars and the 
ISM after birth. 

A naive explanation for the present day systematic behaviour 
of the stellar mass fraction is the accumulation of stars, born in the 
periphery, in the central cluster. If such a process was possible, it 
would boost the observed fracy^^ in the cluster while reducing its 
value outside. However this is an implausible scenario. Consider 
the free-fall time (rg) of a star born in the periphery to be pulled 
into the cluster by gravity: 


^ d3/2 

2 yjmu' 


( 2 ) 


where d is the distance from the cluster, Md is the mass of the 
cluster (gas plus stars) and G is the gravitational constant, rg is the 
shortest time-scale for a star to migrate towards the cluster, ignor¬ 
ing friction, virial velocities, rotational motions and the existence 
of any matter outside of the cluster that will exercise gravitational 
pull. For the observed Md of ~10® M© 

(= Af*<i 5 pcATcat frac^\ = 2300 • 4.3/0.1) and a distance of 
30 pc (the distance between region #1 and #3 in Fig. [i^ , Tft 
~9 Myr, considerably longer than the ages of the PMS stars. 
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Figure 12. Distinct regions within the NGC 346 complex as observed in 
a number of key tracers: CO(J=2-l) integrate intensity map jRubio et al.| 
|2000| panel a), H« emission jSmith, Leiton & Pizarro||2000| panel b), 
young star density map based on |Gouliermis et al.H2006| panel c) and the 
dust surface density (this work, panel d). All maps are shown in linear scale. 
We indicate four regions corresponding to local maxima in CO emission. 


The second effect that could alter the observed fracy^^ over 
time, is gas removal from the central region. When a cluster is 
formed, the gas left over from the star formation process may 
be rapidly removed by stellar winds ([Lada, Margulis & Dearborn | 
|1984[ l. More recent work, however, suggests that the winds (and 
later, supernovae) may escape through low density channels formed 
in a hierarchical ISM (e.g. [Rogers & Pittard|2013[|Dale et al.|2014^ 
causing much of the residual gas to remain. This seems to be the 
case for NGC 346. HST images of NGC 346 in the light of H a 
show the presence of residual gas. Molecular gas also appears to be 
present in the centre of the cluster despite the proximity of O stars 
(Fig. [12^. Chandr a observations reveal a low diffuse X-ray flux 
( [Naze et al.|2002] (, which indicates either that a supernova explo¬ 
sion has yet to occur, or the hot gas has escaped the region. There¬ 
fore gas expulsion cannot have produced the position dependence 
of fracyi^. In further support of this conclusion are the findings of 
|Smith| ( [2012[ l and Smith et al. (in prep). These authors present high 
spectral resolution observations of H a and [O III] to examine the 
ionized gas motions in NGC 346. They find that the gas is quies¬ 
cent. Data over ten slit positions covering the ionized nebula show 
no evidence for velocity splitting, even at the centre of the cluster. 
These results indicate that the gas is not being expelled by stellar 
winds but rather the dominant form of interaction is via stellar ra¬ 
diation to produce a classical expanding H II region. 

An additional argument against the importance of efficient gas 
removal leading to a cavity in the central part of our map comes 
from the Spitzer detected embedded YSO candidates ( [Simon et al.j 
[2007|[^wilo et al.|2013| l. These sources are young (;$1 Myr) and 
several of them are detected towards the central cluster. The num¬ 
ber of detected YSO candidates towards the centre of NGC 346 is 
probably a lower limit, set by the difficulty of detecting individual 
point-sources which exhibit a mid-IR excess in the crowded clus¬ 


ter region against the strong extended emission component due to 
the h ot dust. The SFR d erived from these objects (>3.2 10“^ Mq 
yr“^, Simon et al.|2007| is fully compatible with the average SFR 
for the entire 5 Myr event. If the clearing of the natal gas had taken 
place, one would expect to find few of these very young sources to¬ 
wards the central part of our map and probably also a diminishing 
trend of the star formation rate. 

We conclude that the difference of a factor of ~5 in ob¬ 
served fracyi^, most likely reflects a real difference in the aver¬ 
age SFE in the clustered versus the dispersed star-forming regimes. 
This finding is consistent with theoretical models that assume 
that the star-formation efficiency is proportional to the free-fall 
time, hence denser regions have higher efficiencies (e.g. Elmegreeli| 


[2000[[Krumholz & Tan[2007[[Kruijssen[2012^[Hopkins 2013^. 


7.3 High SFE - linked to low metallicity? 

The SFE of NGC 346 is high even after averaging over ~ 150 pc 
(radius) scales (Fig. |^. As discussed in [Kruijssen & Longmore[ 
( [2014^ , this may be explained by the bias introduced by having 
chosen an active star forming complex as the object of study. These 
authors predict an offset of a factor of ~3 for the size of the regions 
we study , similar to what is observed. Additionally, this high value 
may be related to an above average star-formation efficiency in the 
clustered environment which dominates NGC 346. Nevertheless, it 
is still interesting to discuss the possible role of the lower metallic¬ 
ity for the following reasons, i) [Kennicutt[ P008| l shows that low- 
metallicity galaxies like the SMC, as a whole, lie a bit (factor of 
~3) above the general SK-relation. ii) The largest radius we mea¬ 
sured is 150 pc. It is not obvious that by averaging to a larger dis¬ 
tance, the derived SFE will decrease significantly - as it is also not 
obvious that it will stay high - because the next prominent star¬ 
forming complex (LHA 120-N76) is only about 300 pc projected 
distance away from NGC 346. 

Because the SFE within individual star-forming clouds is a 
critical parameter for theories of galaxy evolution and galactic- 
scale star formation, with observable consequences on for example 
the galaxy luminosity function, there are many recent studies that 
have tried to assess whether the conversion of gas to stars depends 
on local ISM conditions and/or global galaxy parameters. [Sain-[ 
[tonge et ar] ( [201l| l investigated the global gas content and star for¬ 
mation rates in ~200 galaxies in the redshift range between 0.025 
and 0.05, finding a clear decrease in the star formation efficiency 
with increasing stellar mass. Resolved studies of nearby galaxies, 
meanwhile, have shown that the star formation efficiency on ~kpc 
scales varies with several ISM and Galactic-scale parameters, in¬ 
cluding stellar mass, dust-to-gas ratio, rotation velocity, orbital time 
and Hubble type (e.g. [Leroy et al.|20T3| l. 

While the influence of any one parameter is difficult to iso¬ 
late, a general conclusion that can be drawn from these studies 
is that low mass, low metallicity, late-type galaxies appear to be 
more efficient at forming stars from their molecular reservoir than 
high mass, metal-rich, early-type disc-galaxies. This could be in 
part a direct effect of metallicity: [Dib[ j20 11 ^ show that low metal¬ 
licity star forming regions may lie above the SK-relation for nor¬ 
mal metallicity galaxies by a factor of 0.3dex as a consequence of 
higher SFE caused by less effective removal of the gas from the star 
forming region by stellar winds at low metallicity jPib et al.[201 l[ l. 
The observed trend may also be due to an increasing fraction of 
molecular material located outside bound star-forming structures as 
the stellar mass and hence ISM pressure increases, although other 
effects are likely to contribute to the observed trends. 
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Our results in Sect. |5.1| are consistent with the picture emerg¬ 
ing from lower resolution studies that molecular gas in low- 
metallicity, low-mass systems is somewhat more efficient at form¬ 
ing stars. In our case, we find that the mass of young stars in NGC 
346 is high relative to its gas content, even after we average our 
measurements over scales that approach the spatial resolution of 
nearby galaxy surveys. However, we caution the reader to realise 
that the presented analysis of NGC 346 alone does not allow us to 
derive the global SMC star-formation efficiency nor conclude that 
it is above average. 


8 SUMMARY 

We present a detailed comparison of the surface density of young 
(PMS-l-UMS) stars (E*) and dust surface density (Edust) across 
the face of the most prominent star forming complex NGC 346 
(N66) in the SMC. E* is derived from individual stars identified in 
the colour-magnitude diagram based on HST photometry. Edust is 
obtained via radiative transfer modeling of the infrared to submm 
spectral energy distribution. We measure these quantities in 115 
independent pixels of 6x6pc^ out to distances of ~50pc from the 
main cluster, covering an area of >4000 pc^. 

We find a correlation between E* as a function of Edust with a 
considerable scatter. A power law fit to the data yields a steep rela¬ 
tion with an exponent of 2.6±0.2. We find that sight-lines towards 
the central ;$15pc exhibit systematically high values of E*/Edust 
by a factor of ~5 compared to the rest of the complex. 

We convert E* to average star formation rate (SFR) surface 
densities (Esfr in Mq yr^'^pc”^) using 1) the observed mass 
function to derive the total stellar mass formed and 2) the time span 
over which these young stars have formed (5 Myr). This Esfr is 
much less biased to the contribution of high-mass stars than lumi¬ 
nosity based SFR tracers. The derived total SFR (4±M0“® Mq 
yr“^) is consistent (within 10 per cent) with Esfr estimated from 
the extinction corrected H a emission integrated over the entire H a 
nebula. 

On small scales (6x6pc^) the SFR estimated from indirect 
tracers (ionized gas or dust emission) compares very poorly with 
the SFR based on star counts. Across most of the studied region, 
H a systematically underestimates the local SFR, reaching a factor 
of 10 discrepancy in the central part of the complex. This discrep¬ 
ancy is due to ionizing photons escaping the local volume where 
the young stars are counted. 

We convert Edust to gas surface density (Egas) using a con¬ 
stant gas-to-dust mass ratio (rgd = 1250). This value is calibrated 
on our dust SED fitting method to yield the measured gas mass 
(from HI and CO) on a large scale, within a 500 arcsec radius cen¬ 
tred on NGC 346. On this scale most of the gas mass is contained in 
the atomic phase and the influence the uncertain conversion factor 
from CO intensity to H 2 column density is small. 

Using the derived Egas and Esfr, we place the pixels on 
a Schmidt-Kennicutt (SK) diagram. We find that individual pix¬ 
els fall systematically above the SK-relation for integrated disc- 
galaxies by on average a factor of ~7. This offset decreases as 
larger beams are used, consistent with the prediction from |Krui-| 
[jssen & Lon g mor^ (|2014|l. 

The observed offset is reminiscent of the results of IHeiderH 
[man et ^ ( |2010[ l and |Wu et al.| ( [20T0l l who show that the mea¬ 
sured SFR lies above the galaxy averages when ‘zooming-in’ to 
parsec-scale on star forming clumps in the Milky Way. This is 
probably caused by less dense gas, which is inefficient in forming 


stars, that is included in the galaxy-scale averages but not included 
when ‘zooming-in’ on individual clumps. A main difference with 
our analysis is that, even though we analyse 6x6 pc^ sized regions, 
we cover the full star-forming complex and find that the entire ac¬ 
tive region (~50pc) lies consistently above the relation. In other 
words for NGC 346 one needs to ‘zoom-out’ beyond 50 pc to probe 
this less dense, non-star forming gas. The average over a larger area 
(90 pc radius) lies closer to the SK-relation but remains high by a 
factor of ~3. 

We present a map of the observed stellar mass fraction 
(frac^^ = E* (Egas -l-E*)“^). frac^i^ is systematically high 
(~10 per cent) within the central 15 pc of the complex and sys¬ 
tematically lower outside (2 per cent). We interpret this differ¬ 
ence as the results of a higher star formation efficiency in the cen¬ 
tral zone. The inner 15 pc is dominated by young stars belonging 
to a centrally condensed cluster, while the outer parts are dom¬ 
inated by a dispersed population ^Gouliermis, Hony & Klessen] 
|2014[ l. Therefore, the observed behaviour could reflect a change in 
star-formation efficiency between clustered and non-clustered star- 
formation. 
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